Dynamic critical behavior of the XY model in small-world networks 
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The critical behavior of the XY model on small-world network is investigated by means of dynamic 
Monte Carlo simulations. We use the short-time relaxation scheme, i.e., the critical behavior is 
studied from the nonequilibrium relaxation to equilibrium. Static and dynamic critical exponents 
are extracted through the use of the dynamic finite-size scaling analysis. It is concluded that the 
dynamic universality class at the transition is of the mean-field nature. We also confirm numerically 
that the value of dynamic critical exponent is independent of the rewiring probability P for P > 0.03. 
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I. INTRODUCTION 

In recent years, there has been a surge of activity in 
the field of complex networks among statistical and in- 
terdisciplinary physicists pj. Quite naturally, various 
spin models of statistical mechanics have been studied on 
an underlying complex network 0. [1 These studies 
serve a twofold purposes: Firstly, they aid studies of the 
static network structure. In many real-world situations, 
the network structure is an underlying infrastructure for 
a dynamical system, and non-trivial effects can emerge 
from the interplay between the dynamical system and 
the network [5j . Secondly, such studies of spins systems 
on complex networks can illuminate the properties of the 
spin-model itself in certain extreme situations. For exam- 
ple, both the Ising and XY models can display a critical 
behavior similar to high dimensional regular lattices with 
a very low density of couplings (or edges in the network) 
between spins 0, • 

One of the most central complex network models is 
the Watts and Strogatz (WS) model of small-world net- 
works 0. Briefly, this model is controlled by a parameter 




FIG. 1: The construction of a Watts-Strogatz model network. 
One starts from a regular one-dimensional lattice (a). For 
every vertex (we consider the black vertex specifically) one 
goes through the edges on one side (for the black vertex these 
are dashed) . Then, with probability P one detaches the other 
end (b) and re-attach it (c) with the condition that no loops 
(edges starting and ending at the same vertex) or multiple 
edges must be formed. 



P (the "rewiring probability" ) , and by tuning P from 
to 1 one goes from regular to random networks. The 
interesting region is that of intermediate P where the 
network is clustered (has a high density of short circuits, 
or more specifically, triangles) and a logarithmically in- 
creasing average path-length (the path length of a pair 
of vertices is the smallest number of intervening edges). 
In the XY model, each vertex is associated with a two- 
dimensional spin-angle. The XY model has mostly been 
used to study phase transitions in superconductors and 
superfiuids, while it was also applied to e.g., the forma- 
tions of bird flocks Q ■ The static properties of the XY 
model in the WS network have been studied in [j| , where 
critical exponents characteristic of a mean-field transition 
have been found at any nonzero value of P. In the present 
paper, we study the dynamic critical behavior of the XY 
model on the WS small-world network with focus on the 
dynamic critical exponent. 



II. XY MODEL ON WS MODEL NETWORK 

In the WS model for the small- world network 0, a 
regular network is first constructed by arranging iV ver- 
tices in a one-dimensional circular topology and connect- 
ing each vertex to 2k neighbors. Then one goes through 
each edge one at a time, and with the rewiring probabil- 
ity P detaches the far-side of the edge and reconnect it to 
a randomly chosen other vertex (with the restriction that 
loops and multiple edges must not be formed). In this 
manner, a small-world network with the size N is con- 
structed with the model parameters k and P. This pro- 
cedure is illustrated Fig. ^ The former parameter (k) is 
not believed to give any significant change of the network 
structure for k > 1, and thus we fix k = 3 throughout 
the paper. 

The XY model consists of planar spins interacting 
through the Hamiltonian 
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where 8i G [— 7T, n) at vertex i is the spin angle, cor- 
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responding to the phase of the superconducting order 
parameter in the Ginzburg-Landau theory of supercon- 
ductivity. The coupling matrix J^- is given by 



Jji — 



J, if [i, j) is an edge, 
0, otherwise. 



(2) 



For example, in the XY model on a two-dimensional 
square lattice where only nearest vertices interact, we 
have Jij = except when i and j are nearest neighbors. 
For convenience, we measure the temperature T in units 
of J/ks- 



III. SHORT-TIME RELAXATION METHOD 
AND SCALING ANALYSIS 

To investigate the dynamic critical behavior of the XY 
model on the WS network we use the so-called "short- 
time relaxation method" , which utilizes the relaxation 
behavior of the system towards equilibrium from the a 
nonequilibrium initial state. By use of this method sev- 
eral critical exponents have been successfully determined 
for the Ising model |E , for unfrustrated and fully frus- 
trated Josephson Junction array s IT(T| . and for the clas- 
sical Heisenberg spin XY model The major advan- 
tage of the short-time relaxation method (compared to 
dynamical simulations in equilibrium) is the running time 
saved from the avoidance of equilibration. 

The Monte Carlo (MC) scheme is based on the Hamil- 
tonian (0 and the standard Metropolis local update al- 
gorithm |l2t lil^ . The key quantity we measure is [jj 0] 



Q(t) 



N 



sgnf y^cosflj(t) 



(3) 



where the time t is measured in units of one MC sweep, 
( • • • ) is the average over different time sequences from the 
same starting configuration, and average over different 
network configurations, denoted by [•••], should also 
be taken. Here the sign function sgn(cc) measures the 
sign (±1) of x. The initial configuration is chosen as 
9i(0) = 0, giving Q(Q) = 1, and Q(t — > oo) = since in 
equilibrium (t — > oo) there is no preferred angle direction. 
We chose the trial angle 89 = 7r/6; the motivation is that 
it is sufficiently small in order to obtain good convergence 
rate of the quantity we measure while it is big enough to 
make simulations fast 0] ■ 

In order to obtain the dynamic critical exponent and 
detect the phase transition, we use the finite-size scaling 
of the quantity Q. Close to the critical temperature T c 
one expects that in a finite-sized system the characteristic 
time r scales as r ~ N z , while the ratio of the correlation 
volume £y ~ \T — T r \~ v to the system size N gives the 
second argument of the scaling function 0, 0, [3, : 



Q(t,T,N) = F^t/N z 



(T - T C )N 



i/p 



where F(xx,X2) is the scaling function with the property 
F(0, X2) = 1. At T c , where the second scaling variable 
vanishes, the dynamic exponent z is easily determined 
from Eq. (0J by the requirement that the Q(t) curves 
obtained for different sizes of the networks collapse onto 
a single curve when plotted against the scaling variable 
tN~ z . It is also possible to determine T c from Eq. Q) by 
applying an intersection method: Starting from the fully 
phase ordered nonequilibrium state, Q decays from I to 
as time proceeds. For times t where < Q(N, T, t) < 1 
we can fix the parameter a = tN~ z to a constant for 
given N and z. Then Q has only one scaling variable 
(T - T^N 1 ^ and can thus be written as 



Q a (T,N) = F\a, {T~T C )N 



jl/u 



(5) 



If now we plot Q with fixed a as a function of T for 
various N, all curves should have a unique intersection 
point at T — T c . Finally, we can check the consistency 
by using the full scaling form to collapse the data for 
different temperatures and networks sizes onto a single 
scaling curve in the variable (T — T^jN 1 '" at fixed a = 
tN~ z . In addition, this is a consistency check of the value 
of the static exponent v. 

To discuss the finite-size scaling in more detail, the 
form Eq. JSJ is based on the assumption that there is 
only two length scales in the system: the network size 
N (or the number of vertices in the network) and the 
correlation volume £y diverging at T c . However, it is 
known that in the small-world network there is an addi- 
tional spatial length scale related with the distance be- 
tween shortcut endpoints, given byC= (kP)- 1 |I3. Ac- 
cordingly, in the presence of the three competing scales 
(N, £y, and Q, the finite-size scaling function should 
take the form x(t/N z ^ v /N,(/N) [l| Here, we 

aim to use sufficiently large systems with N much larger 
than ( (but, as we will see, this is difficult for small 
P), where x(t/N z , £v/N, C/N) may be approximated as 
x(t/N z , £,v/N, 0). This leads to the above mentioned 
scaling forms JIJ and JSJ without (. 



IV. SIMULATION RESULTS 

We exemplify the critical behavior of Q for WS model 
networks with P = 0.2. This value is quite representative 
for all P values of our simulations, but (as we will discuss 
later) small P requires larger system sizes, longer times 
series, and more averages. Figure shows the finite- 
size scaling of the short-time relaxation given by Eq. Q 
which at T c turns into the simple form: 



Q(t)=F(t/N s ,0) 



(0) 



(4) 



with the only one scaling variable t/N z . In Fig. [2] (as 
well as all other P and T values) we have performed a 
sample average over 100 independent runs for 200 differ- 
ent network realizations. Instead of leaving both z and 
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FIG. 2: Short-time relaxation of Q for P = 0.2 at T = T c = 
2.23. Q is shown as a function of the scaling variable tN~ z . 
z = 0.52(1) is found at the best data collapse (see the Ap- 
pendix) . 



FIG. 4: The dynamic critical exponent z as a function of the 
rewiring probability P. The inset showing T c as a function 
of P is consistent with Fig. 4 in Ref. Q. The dashed line is 
z = 0.54. 
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FIG. 3: Finite-size scaling of the short-time relaxation of Q 
with P = 0.2 and at T = 2.18, 2.20, 2.21, 2.23, 2.25, and 
2.27. The inset an intersection plot with fixed t/L z — a and 
(z,a) — (0.52,3.0); this is consistent with T c m 2.23, while 
the main part of the graph displays the full scaling of Q — 
F(t/N E , (T - T c )L 1/i? ) with the mean-field value of v = 2 Q. 



and z = 0.52. In some cases (typically for small P val- 
ues) the T c has to be slightly altered (from the values of 
Ref. 01) to get both the collapse and intersection plots 
of Fig. correct. We then use z and T c estimated as 
above to make the full scaling plot for Q as displayed in 
the main part of Fig. [3] A very smooth collapse here 
is obtained with v = 2.0 which is again consistent, with 
Ref. 0. 

The procedure described above for P = 0.2 is then 
repeated for various values of P to obtain Fig.0J As one 
can see, except for P < 0.03, z = 0.54(3) [2(j throughout 
the broad range of P. We believe that the nature of 
the transition (and hence z) is independent of P for all 
P > 0. The larger values of z for small P is a result 
of a failure of the assumption that we can neglect the 
length scale £ since N 3> £(~ l/P) cannot be valid for 
small P. The inset of Fig. 0] displays the dependence 
of critical temperature T c (obtained as discussed above) 
upon P and is consistent with what has been obtained 
from static MC 0. 



T c as free parameters, we use T c obtained from static MC 
simulations |3j. Figure[3displays the best collapse onto a 
single curve in a broad range of the scaling variable tN~ z 
with z = 0.52(1) where the number in the parenthesis is 
the error in the last digit (how z is obtained is described 
in detail in the Appendix) [l9j . Just as for static quan- 
tities [3j the obtained z is consistent with higher dimen- 
sional regular lattices (d > 4 to be precise), where 2 = 0.5 
is expected |l6j . 

However, the above method presumes a priori knowl- 
edge of T c . To check out the consistency of determination 
of T c one can use an intersection method described above 
in Sec. IIIII In the inset of Fig. we display Q as a func- 
tion of T for different network sizes N with a fixed value 
of a — tN~ z in the first argument of the scaling form in 
Eq. We find a unique crossing point at T — T c = 2.23 



V. SUMMARY 

In conclusion, we have studied the dynamic critical 
behavior of the XY model on WS model networks by 
means of dynamic Monte Carlo simulations. We have 
used the short-time relaxation method, based on the re- 
laxation from a nonequilibrium state, and determined 
the critical temperature T c , the dynamic critical expo- 
nent z, as well as the static correlation- volume exponent 
v. The dynamic critical exponent was determined to be 
z = 0.54(3) for the networks with rewiring probability 
P > 0.03, while the static critical exponent was found to 
be v k, 2.0. We believe that this result will hold for any 
P > but that the system size needed to confirm this di- 
verges as P — > 0. The exponent as well as two others, 
critical exponents a and (3 of the specific heat and magne- 
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tization respectively have been obtained in Ref. [3j. The 
obtained values v = 2, (3 = 1/2, and a = 0, which also 
have been shown to be independent from the value of P, 
establish the mean-field nature of the transition in XY 
model on WS networks. The result of the present paper 
support this picture and since the upper dimensionality 
for the mean-field theory is d = 4, one can conclude that 
the phase transition in XY model on WS networks is in 
the same universality class as a regular lattice of dimen- 
sionality d > 4. 

An interesting observation is that for a regular hyper- 
cubic lattice this behavior requires a number of edges 
larger than 8 A, whereas for in our simulations we have 
much fewer (3N) edges; and most probably k — 2 (giving 
2N edges) gives the same behavior. We also note that 
there is no additional critical behavior induced by the 
WS model other than the transition from linear ( "large- 
world" ) to logarithmic (small- world) behavior in diame- 
ter as P becomes finite. 
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APPENDIX A: DETERMINATION OF z 

This appendix concerns the estimation of z from data 
collapses as illustrated in Fig. The problem we are 



faced with is that we are looking for a collapse over a 
large range of x — tN~ z , and that the functional form of 
Q(x) is not easily expressed on a closed form or in low 
degree series expansions. To get around this problem we 
partition the x-range in N scg segments Xi, 1 < i < -/V se g, 
and fit a line (a^ + biX, x £ Xi) to the Q point-set within 
each segment (cf. Hi]). Then we sum the square of the 
deviations from the lines: 



A ( z "')= E E (QW-ai-hx) 2 , (Al) 



where Xi(z') is the set of all numerical values of tN~ z 
and thus depends on the value of z' chosen. (Note that t 
and hence x are discrete variables.) Now it is clear that 
if the segmentation can be done so that Q can be rea- 
sonably well approximated by the line segments ai + b{X, 
i.e., if Q(x) is smooth enough, then z = min^/ A(z') will 
converge to the correct value as the number of samples 
and A seg are increased. 

The remaining consideration is how to choose the seg- 
mentation. In general, one needs the segments large 
enough to get a small error in the linear regression, and 
small enough for the line-segment approximation to be 
feasible. In practice the method seems to be rather in- 
sensitive for the partition method. We choose to parti- 
tion the whole range of x in segments of equal length, 
with A seg = 30. The minimization of A is conveniently 
done by a Newton- Raphson method [23. The error in z 
is calculated by jackknife estimation 23]. 
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